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Abstract 

We investigate how simultaneously recorded long-range power-law correlated multi- 
variate signals cross-correlate. To this end we introduce a two-component ARFIMA 
stochastic process and a two-component FIARCH process to generate coupled frac- 
tal signals with long-range power-law correlations which are at the same time long- 
range cross-correlated. We study how the degree of cross-correlations between these 
signals depends on the scaling exponents characterizing the fractal correlations in 
each signal and on the coupling between the signals. Our findings have relevance 
when studying parallel outputs of multiple-component of physical, physiological and 
social systems. 



Many empirical data are characterized by long-range power-law auto-correlations 
as well as by long-range cross- correlations. Such scale-invariant organization 
in both auto-correlations and cross-correlations can be observed either for the 
data variables or their absolute values 



Scale- invariant power-law auto-correlations in stochastic variables can be mod- 
eled by the fractionally autoregressive integrated moving-average process (ARFIMA) 
fMTTl: 



Preprint submitted to Elsevier 



16 December 2008 



oo 

Xt = Yl 0'n{d)Xt-n + et, 
n=l 



(1) 



where d G (—0.5, 0.5) is a scaling parameter, et denotes independent and 
identically distributed (i.i.d.) Gaussian variables with (e*) = and (e^) = 1, 
an{d) are the weights defined by a„((i) = dV{n — d) / {V{1 — d)V{n + l)), where 
r denotes the Gamma function and n is the time scale. We denote the auto- 
correlation function for Xt as A{xtiXt-n) = A{n). For c? = the generated 
variable Xt becomes random. 

To account for power-law cross-correlations between two variables Xt and yt-, 
where each variable is itself power-law auto-correlated, we propose a two- 
component ARFIMA stochastic process defined by two stochastic variables Xt 
and Ut- Each of these variables at any time depends not only on its own past 
values but also on past values of the other variable: 



Xt = [WXt + {l-W)Yt]+et, 
yt = [{l-W)Xt + WYt] + ~et, 

oo 

-^t = X/ '^n{,di)xt-m 
n=l 

oo 

Yt = Y^ anid2)yt-n, 

n=l 



(2a) 
(2b) 

(2c) 
(2d) 



where et and it denote i.i.d. Gaussian variables with (e^) = (q) = and 
(e^) = (e^) = 1, a„((ii) and a„((i2) are the weights defined in Eq. (1) through 
the scaling parameters di and d2 (0 < di^2 < 0.5), and W is a. free parameter 
controlling the coupling strength between xt and yt (0.5 < W < 1). We de- 
note the cross-correlation function between xt and yt as C{xt,yt-n) = C{n). 
For different values of a different degree of cross-correlation between the 
variables Xt and yt is observed. For example, for the case when W = 1, the 
process defined in Eqs. ( 12a|) - (l2d|) reduces to two decoupled ARFIMA processes 
defined in Eq. ([1]). Thus, when W = 1 the long-range cross- correlations be- 
tween Xt and yt vanish, while both xt and yt remain long-range power-law 
auto-correlated. 



In Fig. 1(a) we show segments of the time series Xt and yt generated by the 
process defined in Eq. (2a)-(2d) with parameters W = 0.8 and di = d2 = 0.4. 
Both variables exhibit a very similar comovement. In Fig. 1(b) we show the 
auto-correlation functions A{n) for Xt and yt, as well as the cross-correlation 
function C{xt,yt-n) = C{n). These three curves practically overlap [Fig. 1(b), 
three top curves]. We also show the same correlation functions for = 0.8 
and di = d2 = 0.3 [Fig. 1(b), three bottom curves]. Generally, when the 
coupling parameter W is kept fixed, the stochastic process we introduce in 
Eq. (2) generates stronger cross-correlations for larger values of the scaling 
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Fig. 1. (a) Time series xt and yt for the process defined in Eqs. ([2a|) -(WSigv) where 
W = 0.8 and di = d2 = 0.4. The time series xt is vertically shifted for clarity. 
Both Xt and yt exhibit apparent comovement, indicating a high degree of cross-cor- 
relation, (b) Log-log plots of the auto-correlation functions A{n) for xt and yt, and 
their cross-correlation function C{n) for the two-component ARFIMA process with 
W = 0.8 and di = d2 = 0.4 (top three curves), and with W = 0.8 and di = ^2 = 0.3 
(bottom three curves). For decreasing values of the scaling parameters di and d2 
both the auto-correlations and cross-correlations decrease, leading to smaller values 
of A(n) and C{n). 




Fig. 2. Cross-correlation function C{n) before Fourier phase-randomization proce- 
dure for the time series xt and yt shown in Fig. 1 (open symbols). After Fourier 
phase randomization of xt and yt the cross-correlation function virtually disappears 
(filled symbols) for any value of di and d2- 

parameters di and d2. 

Motivated by the fact that for linear processes the auto-correlation function 
does not change under randomization of the Fourier phase [T3|[T^ , we next test 
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Fig. 3. Cross-correlation function C(n) between time series xt and yt generated by 
the process in Eqs.(2) for varying values of W and di = d2 = 0.4. The cross-correla- 
tion function has highest values for W = 0.5, and tends to zero for W approaching 
1. When W = 1, Xt and yt become two decoupled ARFIMA processes. 

how this phase-randomization procedure affects the degree of cross-correlation 
between Xt and yt- First, we perform a Fourier transform of the original time 
series, e.g. Xt, preserving the Fourier amplitudes but randomizing the Fourier 
phases. Then, we perform an inverse Fourier transform and obtain a surrogate 
(linearized) time series Xt- Applying this phase-randomization procedure to 
both time series Xt and yt generated by the two-component ARFIMA process 
in Eq. (2), we calculate the two auto-correlation functions for Xt and ijt, 
as well as their cross-correlation function C{xt,yt-n)- As expected, the auto- 
correlation functions remain unchanged after Fourier phase randomization, 
but the cross-correlation function C{xt,yt~n) completely vanishes [Fig. 2]. 

Next, we investigate the case when the scaling parameters di and d2 are 
fixed, while the coupling parameter W varies. In Fig. 3, we show how the 
cross-correlation function changes for different values of W and for fixed 
di = d2 = 0.4. The closer the value of the parameter W to 1, the weaker the 
cross-correlations {W = 1 corresponds to the case of two decoupled ARFIMA 
processes). 

Next we analyze how the degree of power-law auto-correlations changes when 
varying the parameters W, di, and d2 in Eqs. fl2aj) - fl2dl) . To quantify the auto- 
correlations we employ the detrended fluctuations analysis (DFA) method. 
We estimate the rms fluctuation function F{n) for different time scales n 
[TS|16|17|18|19] . A power-law dependence of F{n) on the time scale n — 
F{n) oc n", where a is the correlation exponent — indicates presence of power 
law auto-correlations. In Fig. 4, we show the DFA scaling curves obtained for 
Xt and yt generated by the two-component ARFIMA process in Eqs. (l2a|) - pd|) . 
where di = 0.4 and d2 = 0.1, and the coupling parameter W varies. For W = 1 
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Fig. 4. DFA scaling curves for the time series xt and yt generated by the two- 
component ARFIMA process in Eqs. (j2ap - ()2dp . where di = 0.4 and ^2 = 0.1. For 
W = 1, Xt and yt are decoupled and thus not cross-correlated, and xt behaves as 
the ARFIMA process in Eq. ([1]) defined only by the scaling parameter di, while 
yt becomes a separate ARFIMA process defined only by the scaling parameter d2- 
For W 1, the scaling properties of xt depend on both parameters di and d2- 
When W = 0.5, the DFA correlation exponent a for xt becomes equal to the DFA 
correlation exponent for yt. The DFA exponent for \yt\ does not depend on W. 

the processes Xt and yt are decoupled and thus not cross-correlated. In this 
case, Xt behaves as a power-law auto-correlated ARFIMA process controlled 
by only the scaling parameter di, with the DFA correlation exponent equals 
a = 0.5 + di = 0.9. Similarly, yt becomes a separate ARFIMA process (de- 
coupled from Xt) which is controlled only by the scaling parameter d2, where 
a = 0.5 + d2 = 0.6. We find that with decreasing value of W (from 1 to 0.5), 
Xt becomes a mixture of two ARFIMA processes and the DFA correlation ex- 
ponent a gradually decreases towards a = 0.6 corresponding to the yt process, 
controlled by parameter = 0.1. In contrast to Xt, for the process yt the DFA 
correlation exponent a virtually does not change with varying the coupling 
parameter W. 

We next consider a separate stochastic process which generates simultaneously 
two time series with power-law auto-correlated absolute values of their vari- 
ables and long-range cross-correlations between these absolute values. Power- 
law auto-correlations in the absolute values of the stochastic variables can be 
modeled by the Fractionally Integrated ARCH (FIARCH) process [20p2] : 





Xt = otet 



(3a) 



oo 



(Jt = Y^ an{d) 



(3b) 



n=l 
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where denotes an i.i.d. Gaussian variable with (e^) = and (e^) = 1, 
and < d < 1/2 and = {l^tD- The sum of the weights a„((i) satis- 
fies Z]^i r(f-d)r(n+i) ~ yielding {at) = 1. While for the time series Xt 
generated by Eq. (1) the autocorrelation function A{xt,Xt-n) is zero for all 
time scales n, for the absolute values \xt\ the auto-correlation function is 
A{\xt\, \xt-n\) = r(l - d)T{n + d)/{T{d)T{n + 1 - d)), which for n » 1 
converges to the power law A{n) ~ n^^^'^'^. 

To account for power-law cross-correlations between the absolute values of two 
variables, where the absolute values of each variable are simultaneously power- 
law auto-correlated, we have previously introduced [21] a two-component FI- 
ARCH process with scaling parameters di and d2'- 



Xt = [Waxt + {l-W)(Tyt]et 
yt = [il-W)a,t + Wayt]it 

di T{n — di) \xt-n\ 



a. 



xt 



E 



n=l 
oo 



a, 



yt 



E 



r(i-rfi)r(n + i) 

d2T{n-d2) \yt-n\ 



\T{1 - d2)T{n + 1) jjy 



(4a) 
(4b) 

(4c) 
(4d) 



where ej and et are i.i.d. variables with (e() = (e() = and (e^) = {ef) = 1, 
W is the coupling parameter controlling the degree of cross-correlations, and 
/ix = {\xt\) and /iy = {\yt\). 

Note, that each of the variables is controlled by a composite volatility — e.g. 
for Xt the composite volatility is Wiaxt + (1 — Wi)ayt [Eq. (Hal) ] — that is 
a combination of two FIARCH volatilities axt and axt [Eq. (I3b[) ]. Stability 
of the FIARCH process is achieved through the condition (a^) = 1. To retain 
stability for the two-component FIARCH process in Eq. (4), the average values 
of the composite volatilities Waxt + (1 — W)ayt and (1 — W)axt + Wayt in 
Eqs. (I4a l) -f j4bl) should be 1. For W = 1 the process in Eqs. f Haj) -( l4dl) reduces 
to two decoupled FIARCH process as defined in Eqs. ( I3bl) -(l3bl). and thus \xt\ 
and \yt\ are not cross-correlated. 

In Ref. [21] we have analyzed the cross-correlation functions between \xt\ and 
\yt\ for the process defined in Eqs. ( I4al) -( l4dl) for varying values of the param- 
eters W, di, and d2. 

Finally, we analyze how the auto-correlations in the absolute values change 
when varying the parameters W, di, and d2- In Fig. 5, we show the DFA 
scaling curves for di = 0.4 and d2 = 0.1, and for varying W. For W = 1, 
the time series Xt and yt are decoupled and so not cross-correlated. In this 
case, Xt is a FIARCH process controlled only by the scaling parameter di, 
and exhibits long-range power-law auto-correlations characterized by a DFA 
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Fig. 5. DFA scaling curves of the time series \xt\ and \yt\ generated by the two- 
component FIARCH process in Eqs. (|4ap - ()4dp . where di = 0.4 and d2 = 0.1. For 
W = \, \xt\ and \yt\ are decoupled and thus are not cross-correlated. In this case, 
xt becomes a separate FIARCH process as defined in Eqs. (3), and the auto-corre- 
lation properties of xt depend only on the scaling parameter di, while yt is another 
FIARCH process with auto-correlation properties depending only on the parameter 
d2- For W ^ 1, the scaling properties of xt depend on both parameters di and d2- 
When W = 0.5, the DFA correlation exponent a for \xt \ becomes equal to the DFA 
correlation exponent for \yt\. Note that the DFA exponent for \yt\ does not depend 
on W. 

correlation exponent a = 0.5 + di = 0.9. Similarly, yt is another FIARCH 
process controlled only by d2, and characterized by a = 0.5 + ^2 = 0.6. We 
find that with decreasing value of W (from 1 to 0), Xt is controlled by both 
parameters di and d2, and the DFA exponent a gradually decreases towards 
the value a = 0.6. At the same time, the process yt which is controlled only 
by the parameter d2 = 0.1 is also characterized by a = 0.6, regardless of the 
values of W. 

The presented modeling approach and findings may have relevance when quan- 
tifying cross- correlations in simultaneously recorded multivariate time series 
of fractal nature. This problem is pertinent to multiple component physical 
[22|23f24j . physiological, social and financial systems. 
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